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1 Introduction 


A Hamiltonian system with n degrees of freedom is called Liouville (or Liouville-Arnold) integrable if it 
is endowed with n functionally independent integrals of motion in involution (including the Hamiltonian 
itself). Some integrable systems, as the harmonic oscillator or the Kepler problem, admit more constants 
of motion than degrees of freedom; they are called superintegrable. Therefore a Hamiltonian with two 
degrees of freedom is said to be superintegrable if it admits two fundamentals integrals of motion, Ii and 
I 2 , that Poisson commute and a third independent integral I 3 . The additional integral I 3 has vanishing 
Poisson bracket with H but not necessarily with Ii and 12 - 

It is known that if a system is separable (Hamilton-Jacobi separable in the classical case or Schrddinger 
separable in the quantum case) then it is integrable with integrals of motion of at most second order 
in momenta. Thus, if a system admits multiseparability (separability in several different systems of 
coordinates) then it is endowed with “quadratic superintegrability” (superintegrability with linear or 
quadratic integrals of motion). 

Fris et al. studied [1] the two-dimensional Euclidean systems, which admit separability in more than one 
coordinate systems and obtained four families of potentials Vr, r = a,b,c,d, possessing three functionally 
independent integrals of motion (they were mainly interested in the quantum two-dimensional Schrddinger 
equation but the results obtained are also valid at the classical level). Then other authors studied 
similar problems on higher-dimensional Euclidean spaces [2]-[4], on two-dimensional spaces with a pseuo- 
Euclidean metric (Drach potentials) [5]-[8], or on curved spaces [9]-[15] (see [16] for a recent review on 
superintegrability that includes a long list of references). 

For some time the studies on superintegrability were mainly concerned with “quadratic superinte¬ 
grability” but recent studies have proved the existence of certain systems endowed with higher order 
superintegrability, that is, with integrals of motion which are polynomials in the momenta of higher order 
than two. We mention the Calogero-Moser system whose superintegrability is related with a Lax equation 
[17]-[19] (this formalism is not considered in this paper) and three important systems that are separable 
but in only one system of coordinates, the generalized SW system system (caged anisotropic oscillator) 
[20]-[22], the Tremblay-Turbiner-Winternitz (TTW) system [23]-[37], and the Post-Winternitz (PW) 
system [38]-[39]; in these three cases two of the integrals are quadratic but the third one is of higher 
order. 

One important point is that although the number of superintegrable systems can be considered as rather 
limited, they are not however isolated systems but, on the contrary, they frequently appear grouped in 
families; for example, everyone of the above mentioned potentials Vr, r = a,b,c,d, has structure of a 
three-dimensional vector space. Now in this paper we consider the idea of one-parameter deformations 
of a given Hamiltonian, that is, families of Hamiltonians depending of a real parameter k that are 
superintegrable for all the values of k (in the domain of the parameter) and that for k = 0 they reduce 
to superintegrable Hamiltonians previously studied. 

The main objective of this article is twofold. First, study the existence of superintegrable deforma¬ 
tions of the four two-dimensional Euclidean systems (some of the results obtained are related with some 
nonlinear systems studied in [40]-[43]) and, second, study the existence of superintegrable deformations 
of the TTW and the PW systems. In these two cases we prove the superintegrability and we obtain the 
explicit expression of the third integral (that we recall is of higher-order) as the product of powers of two 
particular complex functions (this complex formalism is very similar to the approach presented in [30], 
[37], for the TTW system and in [39] for the PW system) 

The plan of the article is as follows: Sec. 2 is devoted to recall the main characteristics of the 
four two-dimensional potentials whith separability in two different coordinate systems in the Euclidean 
plane. Then Sec. 3 and 4 are concerned with quadratic superintegrability and Sec. 5 with higher order 
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constants of motion. In Sec. 3 we first introduce the idea of deformation (depending of one parameter 
k) of a Hamiltonian and then we study four families of superintegrable families endowed with multiple 
separability. Sec. 5 has two parts. In the first part we study a ^-dependent Hamiltonian related with the 
harmonic oscillator, that can be considered as a deformation of the TTW system, and in the second part 
a K-dependent Hamiltonian related with the Kepler problem, that can be considered as a deformation of 
the PW system. Finally in Sec. 6 we make some final comments. 


2 Superintegrability with quadratic constants of motion in the 
Euclidean plane 


Let us denote hy Vr, r = a, b, c, d, the four two-dimensional potentials whith separability in two different 
coordinate systems in the Euclidean plane. 

The two first potentials, 14 y Vb, are related with the harmonic oscillator. They satisfy the equation 
Vxy = 0 and correspond, therefore, to a direct sum of one-degree of freedom systems. 


(a) The following potential 


V, = il)a^{x^+y^) + 





( 1 ) 


is separable in (i) Cartesian coordinates and (ii) polar coordinates. The three constants of motion, 
lai, Ia2,y laS, are given by 


T 2 I 2 2 I 2^3 

Ia2 — Py + Ct y - 3 “ , 

y2 

Ia3 = {xpy - ypx)'^ + 2k2{-)^ + 2k3{-)^ . 

X y 


(b) The following potential 

Vb = {\)a^{'ix'^+y‘^) + ^ + kzx ( 2 ) 

is separable in (i) Cartesian coordinates and (ii) parabolic coordinates. The three constants of 
motion, /hi, Ib2, y hz, are given by 


hi 

Ib2 

hs 


p^ -I- 4a^x^ + 2k3X , 

2 . 22, 2 k2 

Py4- a y H-^ , 

y y2 

{xpy - ypx)Py - a^xy"^ + ^ 2 (^) 


h 

2y2 ■ 


The two other potentials, 14 y Vc, are related with the Kepler problem. 


(c) The following potential 




ksx 

y'^y^x'^ + 2/2 


( 3 ) 
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is separable in (i) polar coordinates and (ii) parabolic coordinates. The first constant of motion is 
the Hamiltonian itself, that is Id = He, and the other two constants of motion, 1 ^2, and are 
given by 


Ic2 

Ic3 


{xpy -ypxf 


2 k 2 X^ 


2 k^x\/x‘^ 


{xpy - ypx)py + 


kix 

+ j/2 


y^ 

2 k 2 X 

y2 


k3{2x'^ +y'^) 


(d) The following potential 


Vd 


ki , [^a;2 +y'^ + x] + y'^ - x] 

- , - + k 2 - , :-h K 3 - , :- 

a/x^ + a/x^ + j/2 


( 4 ) 


is separable in (i) parabolic coordinates (a, &) and (ii) a second system of parabolic coordinares 
(a, /3) obtained from (o, b) by a rotation. The first constant of motion is the Hamiltonian itself, 
that is Idi = Hd, and the other two constants of motion, Id 2 , and Ids, are given by 


Id2 

Id3 


{xpy - yPx)Py + 
{xpy - ypx)Px - 


kix 

+2/^ 

hy 

yx2“+y2 


y [^x2 + 7/2 

/C2- , :- 

X [ a/x^ + — xl 

ko —^- —- 


+ fca 

+ ^3 


y ^x^ + + X 
X [■\/x2 + + x] 

yx2~+^ 


3 Four new superintegrable families endowed with multiple sep¬ 
arability 


Suppose we are given a Hamiltonian H ; then we can construct a new Hamiltonian H as H = XH where 
A is a certain function defined on the configuration space. This new Hamiltonian represents a new and 
different dynamics; for example if H is defined on an Euclidean space then the new dynamics will be 
nonEuclidean. More particularly, we are interested in multipliers A that preserve certain properties of 
H as integrability or separability. For example if H is defined in the Euclidean plane and is separable 
in Cartesian coordinates (x,y) and A is of the form A = l//i, p = /(x) + g{y), then H is separable 
in Cartesian coordinates as well, and if H is separable in polar coordinates (r, </>) and A is of the form 
A = 1/p, p = /(r) + g{(j))/r'^, then H is also separable in polar coordinates. A more strong condition is 
that A must preserve not just separability but multiple separability; this requirement will strongly restrict 
the form of the multiplier. 

Another important property we wish to introduce is ^at the new Hamiltonian H must be a deformation 
of H. By deformation we mean that A, and therefore H, will depend of a parameter k in such a way that 

(i) The new Hamiltonian H(k) must be a continous function of k (in a certain domain of the parameter). 

(ii) When k —> 0 we have A —> 1 and then the dynamics of the Euclidean Hamiltonian H is recovered. 

In this section we study the separability and the superintegrability of four Hamiltonians Hr{K) obtained 
as deformations of the four Hamiltonians H^, r = a, b, c, d. 
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3.1 Hamiltonian Ha^n) 


Let us consider the Hamiltonian Ha 


= \ (pI+pD 

and denote by the following multiplier 

1 


1 2 (^2 

2 X y 


Xa= - , y,a{x,y) = 1 - KT^ 

Ma 


2 2,2 
r = X + y . 


( 5 ) 


Then the new ^-dependent Hamiltonian Ha^n) defined as Ha{K) = XaHa is given by 

1 


SiM = 5 (f^' 




x"^ + 

1 — KT^ 


k 2 h 

{1 — Kr‘^)x'^ (1 —Kr2)y2 


( 6 ) 


The parameter k can take both positive and negative values. In the /« < 0 case the Hain) dynamics is 
correctly defined for all the values of the variables; nevertheless when k > 0, the Hamiltonian (and the 
associated dynamics) has a singularity at 1 — = 0, so in this case the Ha{K) dynamics is defined in 

the interior of the circle = 1 /k, k > 0, that is the region in which kinetic term is positive definite. 


(i) Cartesian separability 

The Hamilton-Jacobi (H-J) equation takes the form 

1 / 1 ^^^/91T^2 /9IT\2-' 1 '-2 

2 vT^ 


/ 1 \ [row \ 2 row \ 

\1 — Kr"^ ) .\ dx ) ”*"Vi9y/J 


,1 2 +y‘^' 

+ 2“ 


(1 — Kr'^)x‘^ (1 — Kr'^)y‘^ 


= E 


so that if we assume that W can be written as H^ = Wx{x) + Wy{y) then we can perform a separation 
of variables and arrive to the following two one-variable expressions 


rdWx^ 
V dx ) 


+ a^x"^ + 


dx 

r/9iT„\2 2 2 

[[-w) 


2fc3 

y2 


2kEx 


2 KEy 


= K + E, 

= -K + E, 


where K denotes the constant associated to the separability. Everyone of these two expressions 
determine a constant of motion. So the ^-dependent functions Jai and Ja 2 given by 




Jal = Px+a X -\ -^ -I- 2KHaX , 


2 k. 


Ja 2 = Py+a y -\ - - + 2KHay 

y2 

are constants of motion satisfying the following properties 


(z) dJal A dJa2 7^ 0 : {'^al : Ja2} — 0 ? (zzz) Hai^rX) — ^ ('^ol “b Ja2^ ■ 


(ii) Polar separability 

The Hamilton-Jacobi (H-J) equation takes the form 

1 r/aH^\2 1 /aiT\2- 

2m -V 9r / r‘^\ drj) ) - 


I 2 2 I ^2 ^3 

+ Oi r H— -— H— . 2 

^ nr\a^ rh Sin"^ ( 


cos^ 


= (1 — Kr‘^)E 


5 




so that if we assume that W is of the form W = Wr{r) + W^{4)) we can perform a separation of 
variables and rewrite the equation as the sum of two one-variable summands 


,fdWr 


dr 


+ a^r'^ — 2Er^{l — nr^) 


-rdW^\^ 2 k 2 2 fc 3 - 

A d(f> ) cos'^(j> sin^ (f>- 


so that the following function 


is a constant of motion. 


Ja3 — 




2 k 2 
2 ^ 


cos 


2k3 

2 j 


sm 


We summarize the results in the following proposition. 


Proposition 1 The K-dependent Hamiltonian Hain) is H-J separable in Cartesian (x, y) and polar (r, 0) 
coordinates and it is endowed with the following three quadratic constants of motion 


1 


(1 — Kr^) 


//I 2\ 2 I 2 2 I 2 2 I 

[{I — Ky )p^ + KX Py + a X -\- 

■ (^(1 - kx^)pI + Ky^pl + a^y'^ + + 


Jal — 

(1 — K r^) 

Ja3 = {xpy - ypxf + 2 k 2 {^Y + 2k3(^Y 


2 ^ 2(1 — ny"^) ^ 2 k 3 X^ 


2 y 2 

2 k 2 y ^ 2/03(1 — Kx ‘^) 

3,2 y2 


We note that the two first functions satisfy the correct limit when k —)■ 0, that is lim„_i.o Jai = lai^ 
i = 1,2. Concerning the third function it is k- independent and it coincides with the original constant 
da3 ■ 


3.2 Hamiltonian 

Let us consider the Hamiltonian Hi, 

^6=2 (Px + pI) + 2 

and denote by At, the following multiplier 


At, = 


Tb 


yb{x,y) = 1- KX. 


Then the new tc-dependent Hamiltonian Hi,{k) is thus given by 


"‘W =5 (F^) ++ 


'4a;^ + 


k 3 X 


{1 — Kx) y'^ {1 — Kx) 


( 7 ) 


( 8 ) 


(i) Cartesian separability 

The H-J equation takes the form 


( " ) 

\(—Y + i 

/51T\2- 

\1 — KxJ 




+ 


/ +y^\ 

V 1 — KX / 


ksx 


{1 — kx) (1 — kx) 


= E 
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so that if we assume that the function W is of the form W = Wx{x) + Wy{y) we can perform a 
separation of variables and arrive to 


rdWx \ 2 


LV dx 

This means that the following two functions 


\ ^ 9 9 1 r / dWy \ 2 „ „ 

J AlOl X -\~ 2,k^x ‘2hiEx H- —j ex y ~\~ 


2 k-2 


dy 


= 2E. 


Jbi = pI + + 2k3X + 2KxHh , Jb2 = + 


2fc2 

y 2 


are constants of motion satisfying the following properties 

(z) dt /51 A (iJ}j2 7 ^ O 5 \_Jhl ; Jh2\ — O 5 (zii) — 2 ('^^1 Jh'l) ■ 

(ii) Parabolic separability 

If we introduce the following change 

{x, y) —)■ (a,b), X = , y = 2ab , 

then the ^-dependent Hamiltonian becomes 


Hb{K) = - 


1 


2 1 - K (a2 - 

and the H-J equation 


_ ( Pa+Pb \ , _ 

52) Va2 + 52; 1 


r 1 o / a® + 5®' 

ot 


70“ + 6“ \ 

-k(o2-52) [ 2 "" Va2 + 52j 


k 2 a‘^ -b-^- 
02^2 + 3 0252 _ 


1 1 


r/aiT\2 /aTT\2 


2 (a2 + b‘^)[\ da J db J \ ^ [2 

also admits separation of variables 

^(dWa\^ , .2„6 , ^k 2 
A da 

The result is that the following function 


rl 


, 7a® + &®\ k2 a'^-b'^l , 2 ;2^^n 


) +a^a® + ^+2A:3a"‘-2(a^-Ka^)i? + +a^&® + ^-2fc36^-2(5^ + K&"‘)i!l 

/ \ l\ 00 / 0 ^ 


= 0 


■hz = viE a^a^+‘^+2k3a^- 2 {a'^ - Ka‘^)Hb{K), 


pI + a^b'^ + ^ - 2 / 036 ^ - 2 ( 6 ^ + Kb‘^)Hb{K) 


is also a constant of motion. 


The following proposition summarizes these results. 

Proposition 2 The n-dependent Hamiltonian Hi,{k) is H-J separable both in Cartesian coordinates {x, y) 
and in parabolic coordinates {a, h), and it is endowed with the following three quadratic constants of motion 


Jbi — 


1 


— KX) 

j _ 2 I 2 2 I 

Jb 2 — Py A- a y -\- 


[pI + Kxpy + a^{Ax + Ky‘^)x H-+ 2k3X^ , 

2 k2 ^ 


Jb3 = {xpy - yPx)Py - • 


"la-tx)’ + (T^ (- A - + ’f))-h\ 
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3.3 Hamiltonian Hc{n) 


Let us consider the Hamiltonian Hr 




^2 , 


y2i/a;2 + y2 


(9) 


and denote by Ac the following multiplier 


Ac = — , /Xc(x,i/) = 1 - - . 

Me r 


Then the new ^-dependent Hamiltonian Hc{k) is 

= I (^) (P' + ^ + 


k 2 r 


+ 


ksx 


— K {r — K)y‘^ {1 — Kr)y 


( 10 ) 


(i) Polar separability 

The H-J equation takes the form 


-1—I 

2 \r — K/ 


-rdW\^ 1 /9VP\2 
2 \r — kJ . \ dr ) ~^r‘^\d(j)) 


ko r 


fcs cos (j) 


—+ - 

r — K {r — K)r sin^ (p {r — K)r sin^ 


= E 


so that if we assume that W is of the form IT = Wr{r) + W^{(j)) we can perform a separation of 
variables and rewrite the equation as the sum of two one-variable summands 


i-w) +2ti-'-2-'(-'-«)El + [(^) +^ + 


2k- 


COS( 

3„- 2 


sin p sin pi 


= 0 


so that the following function 


2 k 2 


Jc2 — P^+ . 2 j 2^3^ 

^ cm ^ Cl 


cosp 


sin^ p sin^ p 

is a constant of motion. 

(ii) Parabolic separability 

The Hamiltonian He takes the following form in parabolic coordinates 


2 (a2 + &2) _ K if- + + (aP + b- 


{aP Eh^) — K 
so that the H-J equation becomes 
-rdWV^ rdW \2 


(a2 + 1 “^) — K 


2k\ -\- ^2 




+ ^3“ 


-52 


a262 


V\ da J \ db J 
and it leads to 


-b2 


2ki + ^2 ^ 2^2 -^ ^3 ^ 2^2 ~ 2((a^ J- b^) — k)E 




2 a^E 


LV db J 

so that the following two functions 




= -kE + K 
= -kE- K 


2 k 2 2^3 ^ 2 ^ 


>/c3a — Pa + 2fci -I-^-^-2a Hc{k) , 



Jczb = Pfc + 2A:i + ^ + ^ - 2b‘^Hc{K) , 
are two constants of motion representing two different ^-deformations of ^c3 

(?) dJcZa ^ dJcSb 0 5 \dc3a ? Jc3b} — 0 5 lim Jc3a — lini Jc3b — ^c3 ■ 

K—>-0 K—>-0 

Nevertheless, as the only difference between them is a term proportional to Hc{k) (with k as 
coefficient), we consider as a more appropiate third constant the following function 

Jc3 — Pa~^ “I-2-2- 2a^ Hc{k) + KHc{bv) 

= - + 2fci + ^ + ^ - 2b^Hc{K) + kHc{k)) 

that in a more detailed way is as follows 

iaVb-bVa)+'lipl-pl) + 2k, (a^-b^) 

+ ^ ( 2 ( 0 ^ + b'^) - + b^)) ■ 


Jc3 — 


1 

q ,2 _|_ 52 _ ^ 


Proposition 3 The K-dependent Hamiltonian Hc{k) is H-J separable in polar coordinates o^nd 

parabolic coordinates (o, b) and it is endowed with the Hamiltonian as the first constant, that is Jci = 
Hfin), and the following two additional quadratic constants of motion 


Jc2 = {xpy - ypxf + 

1 


2 k 2 X^ 


2 kzX\/ 


Jed. — 


r — K 


r {xpy - ypx)py + 2 {xPx - xPy + 2.ypxPy) + kix 
+ (2r - k) + ^ (2 x'^ + y'^ - Kr) 


The function Jc 2 it is K-independent and it coincides with the original constant Ic 2 (the same situation 
we found in the case (a) with Jad)- Concernig Jed, it is clear that it satishes the limit Jed —t led when 

K —)■ 0. 


3.4 Hamiltonian Hd^n) 

Let us consider the Hamiltonian Hd 


„ _ 1 / 2 , 2 \ki Vr + x y/r-x 
kJd — n [Px + Py) ”1-1" ^2-::-1- ks - 


2 2,2 
r = X + y , 


and denote by Xd the following multiplier 


Ad = — , Pd{x,y) = 1 - - . 
Td r 


Thus, the K-dependent Hamiltonian we will study is 


^ (fi^) (P'+Py) + fi^ + 


ki k 2 \/r + X ksy/r — x 


( 11 ) 


( 12 ) 
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(i) Parabolic separability I 

The Hamiltoniano Hdin) takes the following form when written in parabolic coordinates 

^ (a 2 + 62) _ ^ (Pa + p!) + („2 + ^2) _ ^ (^1 + 

SO that the corresponding H-J Equation 


1 

1 


/5IT\2- 

2 ^ 

-1- 62) — K 

K 5a / ^ 



+ 


1 


(a^ + b^) — K 


(^ki + k 2 a + ksb^ = E 


admits separation of variables and it reduces to 


fdWa^ 

V da ) 


T ki T ‘2k2a — ‘2a^E 


fdWby 

LV db ) 


+ ki + 2^36 — 2 b^E 


= —2kE . 


Thus, the following two functions 


Jd2a = pl + ki+ 2k2a - 2a^H ji{k) , 
Jd2b = Pb E ki 2/036 — 2b^Eldi^K), 


that written with more detail are as follows 
1 


Jd2a — 


Jd2b = 


(a^ + b^) — KL 
1 

(o^ + 6^) — K L 


— 6"^p^ + Kp^ + ki(a"‘ — 6^ + k) — 2 k 2 a{r — k) + 2/030^6 
a^pl — b^p^a ~ ^Pb + ki{a'^ — b^ — n) — 2 k 2 ab^ + 2 / 036 ( 0 ^ — n) 


are two independent constants of motion. Nevertheless as the on^ difference between them is just a 
term proportional to the Hamiltonian, that is Jd 2 b — Jd 2 a = 2 k Eldin), we consider more convenient 
to choose the following function 


Jd2 — Jd2b b^Hdi^K) — Jd2a T biHdi^K) 


that takes the form 

■^^2 = q2 ^2 _ ^ (Pa-Pb)) + ^2_^fc2_^ (^l (a2-62)-/c2a(26^-K)+/C3K2a^-K)) 

as the integral representing the K-dependent version of Id 2 - 
(ii) Parabolic separability II 

We can introduce a new system of parabolic coordinates by rotating the original system 
(a,6) (a,/3) ; a= ^(a + /3), 6= ^(a-/3), 

in such a way that the Hamiltonian Eld^n), when written in this new system, it takes the following 
form 


1 


Hdin) = V 


1 


2 (a2 + / 32 )_ 


( 2 . 2 \ 

\Pa+Pft) 


1 


(a2 + — K L 


fci + ^ 2^(0 + /?) + /C3^(a - /3) 


The H-J equation Hd{dW/da, dWjda, a, P) = E is separable as well, and it leads to the following 
two functions 


Jd3a — Pa E ki + V2k20i + V^k^a — 2a^Hd{K), 
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Jd3b — Pp + fci + y/2k2P — — 20^Hd{K). 

that recovering the coordinates (a, b) appear as follows 


JdSa — 

Jd3b = 


1 


(a^ + — K L 

1 

(o^ + 52) — K 


{apb-bpa){apa-bpb)- -{Pa+Pb)'^ -ki{2ab+K)+k2ia+b){b‘^-ab-K)+k 3 {a+b){a^-ab-n) 


{apb-bpa){apa-bpb)+-{pa-Pby-ki{2ab-K)-k2{a-b){b'^+ab-K)+k3{a-b){a^+ab-K) 


This situation is similar to the previous one, that is, they are independent but the only difference 
between them is k times de Hamiltonian, that is J^^b ~ Jd 3 a = “2 k H^{k); therefore we choose the 
following function as the integral of motion associated to the system {a, /3) 


Jd3 = 


1 


q, 2 _j_ ^2 _ ^ 


{apb-hpa ){apa-bpb)-K PaPb 


1 


q,2 _|_ ^2 _ ^ 


(^—2ki ab—k2 b{a^—b^+K)+k3 a(a^—5^—. 


Proposition 4 The K-dependent Hamiltonian Hdin) is H-J separable in two different systems of parabolic 
coordinates, the origial system {a,b) and the rotated system (a,P), and it is endowed with three quadratic 
constants of motion; the Hamiltonian as the first constant, that is Jdi = Hdin), and the following two 
additional integrals 


Jd2 

Jd3 


1 

r — K 

1 

r — K 


r {xpy - yPx)Pv + ^ {xpl - xpl + 2yp,„py^ +kix 

— k 2 {y^/r — X — kV r + + fcs {y^/r + x — nffr — x^ 

r {xpy - yp,,)p,, + ^ [vpI - ypl - 2xp,„p^ -kiy 

— k 2 {x + K)ffr — X + k 3 {x — K)ffr + x 


3.5 Comments 

The previous paragraphs can be considered as rather technical; so now we comment some of the char¬ 
acteristics of these new families of Hamiltonians and we analyze to what physics systems they seem to 
correspond. 


(i) All these “new Hamiltonians” are of the form of the original Hamiltonian multiplied by a factor. 
How was this factor found ? The answer is that every is obtained by imposing separability in 
two different systems of coordinates and this property determines with the only ambiguity of 
the appropriate position and sign of the parameter k. 

It is clear that the factors A^ introduce a change in the geometry of the configuration space. The 
new deformed Hamiltonians Hr{K) describe dynamics in Riemannian manifolds (a particle moving 
in a curved space); the negative sign in front of k has been chosen in order that the sign of n 
coincides with the sign of the curvature. 

(ii) The expressions of the four Hamiltonians Hrin) can be considered as rather involved mainly because 
the nonlinearity affects to both the kinetic term and the potential. Nevertheless the fundamental 
point is that, in spite of their complex aspect, they are directly related with the two fundamental 
superintegrable systems, that is, harmonic oscillators and Kepler problems. 
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The two first Hamiltonians, Ha{K) and are related with the harmonic oscillator. More 

specifically, they are deformations of the nonlinear oscillators Ha and Hf,. Therefore Ha{K) and 
Hi,(k) describe nonlinear oscillators in non-Euclidean spaces. 

The other two Hamiltonians, Hc{k) and are deformations of the nonlinear Hamiltonians Hf, 

and Hfi and therefore they are related with the Kepler problem. For that reason Hc{k) and Hfi{K) 
describe two different versions of the Kepler problem in two different non-Euclidean spaces. 

The following section studies these questions with more detail. 


4 The harmonic oscillator and the Kepler problem on curved 
spaces 

In differential geometric terms, the four Hamiltonians Hri^n), r = a,b,c,d, describe dynamics on non- 
Euclidean spaces. The two first Hamiltonians, Hain) or Ht,{K), are related with the harmonic oscillator 
and the other two, Hc{k) or Hd^n), with the Kepler problem. Now in this section we analyze these two 
important ^-dependent systems from a geometric perspective. 


4.1 The harmonic oscillator on curved spaces 

The following two equations 


cPx 

d?x 

dt'^ 


dx 


1 — Kx"^ \ dt 


a^x 
1 — 


= 0 , 


KX 


dx 


1 — Kx"^ \ dt 


9 

ax 

(1 — Kx'^Y 


= 0 , 


(13) 


(14) 


represent two one-dimensional nonlinear oscillators that can be considered as arising from the Lagrangians 

1 


Lt — — 


^ 1 x^ , ^ 1 9x9 1 x^ 

^ -S' and L 2 = — kx ) — - - -^ 

2 1 — Kx^ 2 1 — Kx^ 2' 2 1 — Kx^ 


with associated Hamiltonians 


and 


TT 1/1 2/ 2 I 0:“^ 

//, = 2(1-“ )P.+ 2 137^' 


Ho = 


1 Px 


1 


2 1 — Kx"^ 2 1 — Kx"^ 


(15) 

(16) 


That is, the two potentials are the same (see Figure 1) but the kinetic terms are different (these two 
oscillators are studied in [44]- [49]). 

The two-dimensional versions of these Lagrangians are 

"‘M = 

= - ( - 

2 \1 — Kr^ 


2 2 

■r^vi I - 


)-y( 


r 

1 — KT^ 


+ Vy- K{xVy - yVxY 


1 

I- x'^ +y'^ \ 

2 ' 

< 1 — K / 
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and 


L2{k) 


-(l-nr^){vl + vl) 


I' \ 

2 V1 — K / 


It is known that a symmetric bilinear form in the velocities (vx,Vy) can be considered as associated to 
a two-dimensional metric ds"^ in IR^. In this particular case, the kinetic term in the Lagrangian Li{k) 
considered as a bilinear form determines the following ^-dependent metric 


ds 


2 


(l — 


(I — dx^ + (1 


Kx^) dy^ + 2k xy dx dy 


The second Lagrangian L 2 {k) represents the harmonic oscillator in a two-dimensional space with a metric 
conformally flat 

ds^ = (1 — Kr‘^){dx^ + dy^). 


The two-dimensional versions of these Hamiltonian are 


and 


Hi(.n) = i ((I -Kr^){pl+pl) + kJ^) + ^ 

2(137^)+ 2" il ^)’ 


where J denotes the angular momentum. They satisfy the correct Euclidean limit 
lim Hi(k) = lim iL 2 (K) = ^ (p^ + Py) + ^ a^x^ + y^) , 


K—^0 


and represent two different harmonic oscillators in two-dimensional curved spaces. 


(17) 


(18) 


(i) The Hamiltonian Hi{k), that has been studied in Refs. [48, 50, 51, 52] (although in some cases 
with a trigonometric-hyperbolic notation), represents the harmonic oscillator in a space of constant 
curvature k (sphere 3“^ with k > 0, and Hiperbolic plane with k < 0). We note that the kinetic 
term includes not only the factor (1 — kt^) but also a contribution of the angular momentum J 
with the curvature k as coefficient. 

(ii) The second Hamiltonian H 2 {k), that is just Hain) with ^2 = ^3 = 0, represents an harmonic 
oscillator in a curved space (two-dimensional space of nonconstant curvature) [40]-[43]. In two 
dimensions the tensor Rabcd only has one independent component which can be taken i?i 2 i 2 

7 ?i 212 = 2 {p2dig2i — 02911 + i 9 i92512 — 9^522) — 9 ef {^11^22 ~ ri2r2i) • 


The result is 


7?1212 — 


2 K 

1 — K 


and the scalar Gaussian curvature is given by 


7?i212 _ ‘2 k 

det[5] (1 — 
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4.2 The Kepler problem on curved spaces 

The following two ^-dependent Hamiltonians 


HKi{n) = ^ {{l-Kr‘^){pl+p^y) + Kj‘^^ -g( ^^ r^ ~) 

and 

represent two different versions of the Kepler problem on curved spaces. They satisfy the correct Euclidean 
limit 


lim Hj^i(k) = lim 
K-J-O K-i-O 


1 

2 


{P'. 


2 I 2\ 

Pv) 


9 

r 


(i) The Hamiltonian Hki{k), that has been studied in Refs. [53, 54], represents the Kepler problem 

in the three spaces of constant curvature k; that is, sphere (k > 0), Euclidean plane (k = 0), 
and Hiperbolic plane (k < 0). Also in this case, when k > 0, the dynamics is restricted to the 
interior of the bounded region < 1/k. 

(ii) The Hamiltonian Hk 2 {k) is just the Hamiltonian of the Kepler problem in Hc{k) or Hdin). In this 
case if k is negative then the Hamilltonian is well defined for all the values k < 0, but when k > 0 
then the dynamics is only well defined in the region r > k in such a way that when k —)■ 0 then the 
dynamics is defined in the whole space. This system is endowed with the following two integrals of 
motion 


J 2 

J 3 


1 


r — K 

1 


r — K 


r {xpy 
r {xpy 


yPx)Py - gx 

ypx)px Egy 


2{r'"-n) ^ ‘^VPxPy) 


that represent the K-dependent version of the (two-dimensional) Runge-Lenz vector. In fact, it can 
be verified that their Poisson bracket is given by 


{J2 , J3} = ^{xpy - ypx) Hk 2 {k) 


that is the same relation characterizing the Runge-Lenz vector in the Euclidean case. 

The metric ds"^, that is also conformal, is given by 

dsl = (1 - -H V). 

We note the factor {1 — k/t) shows a certain similarity with the coefficient (related with the singularity) 
in the Schwarzschild metric. The curvature tensor R 1212 takes the value 

= 2r2(r-«) 

and the Gaussian curvature (scalar) is given by 

_ R 12 I 2 _ K 

det[(7] 2{r—K)^ 
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5 Superintegrabilty with higher orden constants of motion 


We have proved that the two Hamiltonians Ha{K) (related to the harmonic oscillator) and Hc{k) (related 
to the Kepler problem) are superintegrable with quadratic constants of motion. Now, in this section, we 
will prove that they can be considered as particular cases of a more general situation; that is, they admit 
superintegrable generalizations that are separable but in only one system of coordinates (two quadratic 
constants) and are endowed with an additional constant of higher order. 

At this point (as a previous comment to the next two subsections) we notice that a Hamiltonian of the 
form 


1 

_ 1 

( ^ :PI^ 



1 Bi^) 

2 


2 \1 — J 

' + 21 



where A{r) and B{<j)) are arbitrary functions, is H-J separable in polar coordinates (r, (j)) and it is endowed 
with the following constant of motion 

J = pl + B{(j)) = - [r'^pl + A{r) - - Kr‘^)H{K)] . 

This property is true for all the values of the parameter n. 


5.1 A K-dependent Hamiltonian related with the harmonic oscillator and the 
TTW system 


The following potential 


Vuw{r,(j)) — 2 ^ f + . 2^ .n) > 

^ 2r^ Vcos^(m(/)) s\n (m(j))' 


was studied by Tremblay, Turbiner, and Winternitz [23]-[24] and also by other authors [25]-[37]. In 
the general m ^ 1 case it is only separable in polar coordinates (r, (j)) (m must be an integer or rational 
number) and therefore the third integral is not quadratic but a polynomial of higher order in the momenta 
(the degree of the polynomial depends of the value of m). 


Let us consider the K-dependent Hamiltonian 


2 ( 


1 — nr' 


2 

■) {pI + ^), C/am(r,(/)) = 


\ 1 Fmjfl)) 

1 —2(1 —Atr2)r2 


where Fm{(j)) denotes the following angular function 


FM = 


cos^ {mcj)) sin^(TO())) 


and ka and kb are arbitrary constants. It represents a generalization of the Hamiltonian F[a{K) in the 
sense that if m = 1 then we have Flaiin) = Hain). It is clear that this more general Hamiltonian is 
separable but only in polar (r, (j>) coordinates. Therefore it is integrable with the Hamiltonian itself as 
the first integral and a second quadratic constant of motion J 2 associated to the Liouville integrability 


A — nr' 




F^{(l)) 

(1 — nr‘^)r‘^ 


— 7*0 + Pm{4>) ■ 
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Let us denote by Ar and be the complex functions Ar = Ari + i Ar 2 and = B^i + i B ^2 with 
real and imaginary parts, A^a and B^a, a = 1 , 2 , be defined as 


Arl — Pr \/ J 2 : 


Ar2 — 


1 


(1 — r- 


- (1 - 2Kr'^)J-^ = ^ '^2 , 


and 


= J 2 cos{2m(j)) + {kb - ka), B ^2 = P 4 > '/J 2 sm{2m4>). 


First, let us comment that the moduli of these two complex functions, that are constant of motion of 
fourth order in the momenta, are given by 


Ar 


— Ji — 4:kJiJ2 — 4o; J 2 


Bfj) I — J 2 — 2.{ka + kb)J2 + {kb — ka) 


The Poisson bracket of the function Ari with Ham{K) [time-derivative] is proportional to Ar 2 and the 
time-derivative of the Ar 2 is proportional to Ari but with the opposite sign 

{ Arl , Bara (^) } — ^ Ar2 , I Ar2 , Bam (^) } — 2 Arl 5 

and this property is also true for the functions B^i and B ^2 

I B(pi , Bam (^) } — ^ ^ A,i Bfp2 , I Bfp2 , Bam (^) } — ^ ^ A/^ B^i , 

where the common factor A„ takes the value 

Ak = jz - \/~J2 , Ao = ^ \/j 2 ■ 

(1-Kr2)r2 

Therefore, the time-evolution of the complex functions Ar and B^ is given by 

— A-p = z 2 Xf^ Aj -, — B(h = z 2 7Ti . 

at at 


Thus if we denote by Km the complex function B^ then we have 

(a™ s;) = Ar b; + A^ 4* 

= 4 ™”^^ (m{i2X^Ar)B*^ + Ar{-i2mX^B*^)^ = 0 . 


d d 

dt dt 


We summarize this result in the following proposition. 


Proposition 5 The K-dependent Bamiltonian Bam{i^) 


BaM 2 (l-AcrO 


pI 


^ 2( 1 Bm{</>) 

2 — 2(1 — Kr^)r^ 


(23) 


representing a generalization of the Hamiltonian Ha{K), and also a K-deformation of the Euclidean TTW 
Hamiltonian Httw, is superintegrable with two quadratic constants of motion and a third constant of 
motion of higher order. 
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(i) is separable in polar coordinates {r,4>) and it possesses therefore two quadratic constants of 

motion associated to the Liouville integrability 


- (r^^) {Pr + ^) 

•h = p'^+Fm{<P)- 




.2 > 


fiij Ham(.F) admits a complex constant of motion Km defined as 


Km=Af^B;. 


The function Km can be written as Km = J 3 + iJi with J 3 and J 4 real constants of motion. One of them 
can be chosen as the third fundamental integral of motion. 


5.2 A K-dependent Hamiltonian related with the Kepler problen and the PW 
system 


Let us first note that the potential Vc, that in polar coordinates becomes 


9 1/ ^2 fca cos(/) \ 

r Vsin^ (f sin^ <j) / ’ 


can also be written as follows 


t"c = - 


9 

r 


1 / a 
r^ Vcos2((/)/2) 


+ 


/3 


sin"^((/)/ 2 )/ 


k 2 = 2 (q; + ft), ks = 2(/3 - a). 


Therefore, the angular-dependent functions in the potentials 14 and I 4 appear as two particular cases, 
TO = 1 and TO = 1/2, of the general function Fm{4>). It seems therefore natural to conjecture that any 
integrable generalization (or deformation) of the potential 14 must determine a similar generalization (or 
deformation) of the potential Vc. 

In fact there is another interesting system rather similar to the TTW system but that is related, not 
with the harmonic oscillator, but with the Kepler problem (hydrogen atom in the quantum case) 


Vpn,(r, (f) 


-g + lf ^ 

r Vcos^(to(/)) 


_^v 


(24) 


The first study of the superintegrability of this new potential was presented by Post and Winternitz [38] 
by relating Vpw with Vuw making use of the so-called coupling constant metamorphosis transformation 
(Stackel transform) [55]. 

Let us now consider the /t-dependent Kepler-related potential 


He 


2 


lir, </>) 


9 ^ Fmif) 

r — K {r — K)r ' 


(25) 


where Fm{4>) is the same angular function as in the oscillator Ham{F) (and also in this case ka and kb are 
arbitrary constants). It represents a generalization of the Hamiltonian Hc{k) in the sense that if to = 1/2 
then we have LIcm(^)|m=i /2 = Hc{k). It is clear that this more general Hamiltonian is separable but 
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only in polar (r, (j)) coordinates. Therefore it is integrable with the Hamiltonian itself as the first integral 
and a second quadratic constant of motion J 2 associated to the Liouville integrability 


Ji{k) 

J 2 


\r — K/ \ / r — K (r — K)r 

+ 2Fm{(j>) ■ 


We will prove the superintegrability of Hcm{K) using as an approach the existence of a complex factor¬ 
ization for the additional constant of motion. The method, that is in fact a deformation of the formalism 
introduced in [39] for the superintegrability of the PW system, is similar to the one presented in the 
previous section for the Hamiltonian Hamin) (introducing the appropriate changes). 

Now we denote by Ar and be the complex functions Ar = Ari + i Ar 2 and = B^i + i B ^2 with 
real and imaginary parts, Ara and B^a, 0 , = 1,2, be defined as 

Arl =Pr'/^, Ar2 = (^7“) (^-9 r + J 2 + '^{^pI “ y)) . 

and 

Bfjti = J 2 cos{2m<j)) + 2{kb - ka ), B ^2 = P<p sin{2m(j)). 

First, we note that the moduli of these two complex functions are integrals of motion of fourth order 
in the momenta given by 

I Hr P = JiJ2 + {^) Ji—i^gJ2+9^ 

I B(f, P = df — 4(fca + kb)J2 + Mka ~ kb)^ 


The time-derivative [Poisson bracket with iJcm(K)] of the function A^i is proportional to Ar 2 and the 
time-derivative of the Aj .2 is proportional to but with the opposite sign 

{ Hrl , (^) } — Hr2 , { Hr2 , (^) } — Ar\ , 

and this property is also true for the angular functions 

B^^-y , -Hrm (^) } 2 771 B^2 , "{ B(p2 j ^ B(l)\ , 

where the common factor takes the value 

Ak = 7 - 7 — \/~J2 , Ao = ^ \/^ • 

(r-K)r 

Therefore, the time-evolution of the complex functions Aj. and is given by 

—— Hr = i X,^Aj -, -— B{p = i2 ui X^^B^p . 

at at 

Thus, if we denote by Km the complex function Km = Hr^™^ B^, then we have 


dt ™ dt 


h;) = Hr h; + H(2™) Bp 


= 42 — 1 ) (2m(zA„Hr)H; + Hr(-f2mA,5;)) =0. 

We summarize this result in the following proposition. 
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Proposition 6 The K-dependent Hamiltonian 


Hr 




9 


FM = 


r — K Ir — K)r 

h 


(26) 


cos^ {m(f>) sm"^ {m(j)) 


representing a generalization of the Hamiltonian Hc{k), and also a K-deformation of the Euclidean PW 
Hamiltonian Hp^, is superintegrahle with two quadratic constants of motion and a third constant of 
motion of higher order. 

(i) Hcm{F) is separable in polar coordinates {r,(j)) and it possesses therefore two quadratic constants of 
motion associated to the Liouville integrability 


Ji{k) = ( 


r — n 
2 




25 2F„((/)) 


r — K {r — K)r ’ 


J 2 — + 2Fm(<()). 


(ii) Hcm{F) admits a complex constant of motion Km defined as 


6 Final comments 


We observed in the introduction that although the number of superintegrahle systems can be considered 
as rather limited, they are not however isolated systems but, on the contrary, they frequent^ appear 
as grouped in families. Now we have proved the existence of four families of Hamiltonians Hrin), r = 
a,b,c,d, associated to previously known super-separable Hamiltonians Hr- The important point is that 
the multipler A (that is a function dehned in the conhguration space) is a continous function of a parameter 
K so that the superintegrability is preserved for all the values of k and the integrals of motion (and 
therefore the associated symmetries) depend in a smooth way of the parameter. The fact that this 
continous deformations also lead to generalizations of the TTW and the PW systems is certainly a very 
remarkable property. 

We conclude with the following two comments. First, the K-dependent constants of motion are con¬ 
sequence of the existence of K-dependent symmetries; so it would be convenient to study the properties 
of these symmetries from a geometric approach (that is, symplectic formalism and Lie algebra of vector 
fields). Second, it is also convenient to study the quantum versions of these these systems. It is clear that 
these Hamiltonians are systems with a position dependent mass (PDM) and therefore the quantization of 
these systems is not an easy matter. These two points are interesting questions deserving to be studied. 
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Figure 1: Plot of the Potential V = (l/2)a^(a;^/(l — kx^)) as a function of x, for a = 1 and k = 0 (dash 
line), K > 0 (inner upper curves), and k < 0 (lower curves in the outside). 
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